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We study the strong coupling between photons and atoms that can be achieved in an optical 
nanofiber geometry when the interaction is dispersive. While the Purcell enhancement factor for 
spontaneous emission into the guided mode does not reach the strong-coupling regime for individual 
atoms, one can obtain high cooperativity for ensembles of a few thousand atoms due to the tight 
confinement of the guided modes and constructive interference over the entire chain of trapped 
atoms. We calculate the dyadic Green’s function, which determines the scattering of light by atoms 
in the presence of the fiber, and thus the phase shift and polarization rotation induced on the guided 
light by the trapped atoms. The Green’s function is related to a full Heisenberg-Langevin treatment 
of the dispersive response of the quantized field to tensor polarizable atoms. We apply our formalism 
to quantum nondemolition (QND) measurement of the atoms via polarimetry. We study shot-noise- 
limited detection of atom number for atoms in a completely mixed spin state and the squeezing of 
projection noise for atoms in clock states. Gompared with squeezing of atomic ensembles in free 
space, we capitalize on unique features that arise in the nanofiber geometry including anisotropy of 
both the intensity and polarization of the guided modes. We use a first principles stochastic master 
equation to model the squeezing as function of time in the presence of decoherence due to optical 
pumping. We find a peak metrological squeezing of ~ 5 dB is achievable with current technology 
for ~ 2500 atoms trapped 180 nm from the surface of a nanofiber with radius a = 225 nm. 


I. INTRODUCTION 

Strong coupling between atoms and photons is at the 
heart of many quantum information processing proto¬ 
cols including efficient generation of remote entanglement 
[1, 2], quantum data storage and retrieval [3], and QND 
measurements [4]. From a general perspective, strong 
coupling arises when atoms radiate predominantly into 
the electromagnetic field mode that defines the quan¬ 
tum atom-light interface. For an individual atom, the 
strong coupling regime is attained via the Purcell effect, 
whereby the boundary conditions of nearby dielectrics 
and/or conductors enhance radiation into a desired mode 
relative to all other modes. This can be achieved 
with Fabry-Perot cavities (cavity QED) [5] and/or via 
nanophotonic structures engineered such that the radia¬ 
tion is predominantly into a specified mode [6-8]. The 
Purcell enhancement factors for emission into a guided 
or cavity mode scale respectively as Fio/Pvac ~ cto/A 
and Fcav/rvac ~ QX^/V ~ Ta^jA. Here Fvac is the free 
space spontaneous emission rate, uq oc is the resonant 
absorption cross section, Q, V, and IF, are the cavity 
quality factor, volume, and finesse respectively, and A 
is the effective area of the cavity or guided mode that 
couples to the atom. The strongest coupling occurs on 
resonance, and thus much effort has been devoted to de¬ 
veloping the largest possible Fcav and Pid through ultra- 
high-Q, small-volume resonators [-5, 9, 10] and through 
nanophotonic plasmonic [11, 12], metamaterial [13], and 
dielectric [8, 14] waveguides. 

In free space, where there is no Purcell enhancement, 
strong coupling can be achieved via the cooperativity of 
atomic ensembles. This is most naturally implemented 


in a dispersive regime, off resonance, where light elas¬ 
tically scattered from the ensemble constructively inter¬ 
feres to match the mode of an exciting paraxial probe 
[15]. The cooperativity per atom in a typical paraxial 
beam is small, FiD/Pvac ~ cr^/A ^ 10“®. The total 
cooperativity, however, can be significant for sufficiently 
large ensembles, e.g. Na ~ 10^ atoms. The key parame¬ 
ter that characterizes cooperativity is the total resonant 
optical density of the ensemble, OD = NA{cro/A). Such 
strong cooperativity in free space has been employed 
in a variety of applications including quantum memory 
for storage of photonic states [16] and the generation of 
squeezed states of the collective spin of the ensemble via 
quantum nondemolition (QND) measurement [17-20]. 

A particular system that combines the elements above 
consists of cold atoms trapped in the evanescent field of 
the guided mode of a tapered optical nanofiber with a 
subwavelength diameter [21-25] (see Fig. 1). The typi¬ 
cal resonant OD per atom, or OD /Na, in the nanofiber 
{(To/A ~ 10“^) is boosted by orders of magnitude over 
free space for paraxial beams. However, one cannot reach 
the strong coupling regime where Fid is on the order 
of Fvac as is possible in engineered nanophotonic waveg¬ 
uides, such as those arising in photonic crystals [8], where 
atoms can be trapped at positions of peak intensity of the 
field. One can, however, achieve strong cooperativity in 
the dispersive regime with a moderately sized ensemble. 
When compared to free space, all light scattered into the 
guided mode is automatically mode matched, and thus, 
given the relatively large ratio (Tq /A, one can achieve high 
OD with only a few thousand atoms (see Fig. 1). Such 
strong cooperativity opens the door to new regimes to 
create non-Gaussian quantum states of the ensemble [26] 
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FIG. 1. Cooperativity and mode-matching for various atom- 
light geometries, (a) The beam area at the waist of a tightly 
focused beam is closely matched with the atomic scattering 
cross section, but the scattered light of a single atom Is poorly 
mode-matched with the probe, (b) A paraxial beam probing 
a rarefied atomic cloud whose scattered radiation interferes 
constructively in the forward direction, (c) Atoms trapped in 
a ID optical lattices near the surface of an optical nanohber 
interacting with a fiber-guided probe. The tight confinement 
and automatic mode matching that accompanies scattering 
into the guided mode leads to strong cooperativity in the 
atom-photon interaction. 


and potentially to implement nonlinear optics at the level 
of a few photons [27-29]. 

One-dimensional optical lattices in nanofibers based 
on multiple co- and counter-propagating trapping beams 
have been loaded with up to several thousand alkali 
atoms [21, 22]. This has proved a fruitful platform for 
quantum information processing. The anisotropic nature 
of the strong atom-light coupling has been exploited for 
control of internal atomic states [30], enhanced coupling 
into a preferred propagation direction [31, 32], and op¬ 
tical switching [29]. Off resonance, dispersive coupling 
has allowed for non-destructive atom counting [33, 34] 
and storage of fiber-guided light [35, 36]. Recent demon¬ 
strations of photonic crystal cavities fabricated on the 
nanohber [37—39] promise further enhanced atom-light 
coupling. 

In this paper we study the quantum atom-light in¬ 
terface in the dispersive regime for an optical nanohber 
geometry. We focus here on the coupling between the 
atomic spin and light polarization induced by the elastic 
scattering of photons by tensor-polarizable cesium atoms 
trapped near the surface of the nanohber. This provides 
an entangling interaction that can be employed to gener¬ 
ate spin squeezing via QND measurement. Our analysis 
unihes a variety of different approaches found in the liter¬ 
ature, including direct calculation of the dyadic Green’s 
function for photon scattering [6, 11, 40-45] and the 
input-output formalism studied for one-dimensional held 
theories based on Heisenberg-Langevin equations [46-51]. 

The remainder of this article is organized as follows. 
In Sec. II we solve for the mode decomposition of the 
dyadic Green’s function which determines the electric 
held scattered by a point dipole near the surface of the 
nanohber. This allows us to calculate the phase shift 


and polarization transformation for hber-guided photons 
induced by tensor-polarizable atoms in the dispersive 
regime. We connect this with a fully quantum mechani¬ 
cal treatment based on a Heisenberg-Langevin picture in 
Sec. III. The formalism we develop is used in Sec. IV to 
study QND measurement of atoms based on polarization 
spectroscopy. We consider shot-noise-limited atom detec¬ 
tion as well as measurement-backaction-induced squeez¬ 
ing of spin projection noise. We study squeezing of the 
collective pseudospin associated with ensembles of atoms 
in the atomic clock state and calculate its dynamics based 
on a hrst principles stochastic master equation that in¬ 
cludes both the effects of QND measurement as well as 
decoherence due to optical pumping. We conclude with 
a summary and outlook for future research in Sec. V. 


II. DYADIC GREEN’S FUNCTION AND 
INPUT-OUTPUT FIELD RESPONSE 

Given a point particle with tensor polarizability a at 
position r' near the surface of a nanohber, the held at 
frequency luq is given by the solution to the wave equa¬ 
tion, 

[-Vx Vx +n\r)kl] E(r) = -47rfc2d(3)(r_r') a-E(r), (1) 

where fco = wq/c and n(r) is the spatially varying in¬ 
dex of refraction that describes the hber; Gaussian-cgs 
units are used throughout. For an asymptotic input held 
Ein(r), the scattering solution to Eq. (1) is given by the 
Lippmann-Schwinger equation [44], 

Eout(r) = Ei„(r) -h 0^+^(r,r'; wq) • a ■ Eout(r') (2a) 
« Ei„(r) + G(+) (r, r';oJo)-a- Ei„(r'), (2b) 

where in Eq. (2b) we have made the hrst Born approx¬ 
imation valid for weak scattering. The fundamental ob¬ 
ject that fully characterizes the scattered radiation as 
well as the energy level shift and modihed decay rate 
of a scatterer near the dielectric is the dyadic Green’s 
function, G(r, r';a;o). This determines the scattered held 
from a point dipole at r', Escat(r) = G^+Rr, r'; wq) • d, 
and satishes the equation of motion, 

[-VxVx W{r)kl]G{ry-,ujo) = -ATTklS^^\r-r')I, (3) 

where I is the unit tensor. ^ 

The solution for the Green’s function G(r,r';a;o), fol¬ 
lowing from Maxwell’s equations, has been studied pre¬ 
viously [40, 42, 44]. As we are interested here in the 
forward-scattered components that lead to phase shifts 
and polarization transformations, we directly calculate 
G(r, r'; wq) through a decomposition into normal modes. 
A complete set of eigenmodes in the presence of lossless, 
spatially inhomogeneous dielectric are dehned according 
to the procedure of Glauber and Lewenstein [52] . We seek 
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the eigenmodes fr)(r), indexed by rj, that satisfy the ho¬ 
mogeneous wave equation in the absence of sources, i.e., 
Eq. (1) for S = 0 with the eigen-wavenumber fcg —>■ kj^. 
To do so, one defines functions g,,(r) = n(r)f^(r) that 
form a complete basis, as they are eigenfunctions of the 
Hermitian operator, H{ko) = -;^V x V x + kg, 
according to 'H[kg)^^{r) = A,,g,,(r). The eigenvalue, 
Xyj — (bjg — 10 ^ I determines the wavenumber for a 
given mode at frequency uir^ We are interested specif¬ 
ically in the generalized transverse functions satisfying 
V • [n(r)g^(r)] = 0 with eigenvalues A„ 7 ^ 0 [44]. These 
fall into two categories, guided [r] = and unguided 
(j] = v) modes, which together form a complete, or¬ 
thonormal set for transverse vector functions, 


where w(/3) is the frequency of the guided HEn for a 
given (3. 

For z > z' {z < z'), the contribution of the guided 
modes to the retarded (causal) Green’s function is found 
by the usual displacement of the pole on the positive 
(negative) /3-axis into the upper (lower) half of the com¬ 
plex plane. The result for z ^ z' is [ 6 ] 


G^+)(r,r';a;o) 

=27rfy~kes|;3=b/3o 

b,p 


-2 Wo ■' 




Ub/3o,p(rL)Uj^Qp (rDe 


/\JbPo(z-z') 



(9) 


Jd^rg*{r)-gr,'{r) =Jd^rn^{r){*{r)-ir,'{r) (4) 

=XlsM(r)gp(i'') > (5) 

T) p V 


where — v') is the delta function for generalized 

transverse vector fields [53]. It follows that the gener¬ 
alized transverse dyadic Green’s function can be decom¬ 
posed in terms of the eigenfunctions [40, 42] 


G(^)(r,r';a;o) = -47r^ 


1 


Wq - W 2 


( 6 ) 


where 6 = ± indicates the propagation direction, Vg = 
\du!/dl3\p^l}g is the group velocity at wq, and Q{b{z — 
z')) is a Heaviside function enforcing causality for the 
forward- and backward-scattered fields. In the second 
line, we have suppressed the label Pq as it is implicit in 
the definition of the guided modes at frequency uig. 

Radiative properties of a scatterer (the decay rate and 
energy level shift) are determined by evaluation of the 
dyadic Green’s function at the source point r = r' [45]. 
However, for z = z' we cannot close the contour. Instead, 
we expand the resonant denominator in Eq. ( 8 ) with the 
poles moved to yield the retarded (causal) response. 


where the eigenvalues appear as = c^k'^. The sum in¬ 
cludes both guided and unguided contributions. We fo¬ 
cus here on the guided-mode contribution to the Green’s 
function. 

We treat an optical nanofiber of radius a with step- 
index profile. 


n{rj_) = I 


ni 

n2 


r < a 
r > a' 


(7) 


for a silica core (ni = 1.4469) [54] and infinite vacuum 
cladding (712 = !)• For a cylindrically symmetric dielec¬ 
tric the guided modes are f^(r) = Up(r_L)e*^^/v^, with 
indices n = {j, P,p} for the guided mode with propa¬ 
gation constant /3 at frequency ujg = w(/3) and polariza¬ 
tion p. The transverse mode functions are normalized ac¬ 
cording to / d2rj_n2(r_L)u*(rj_)-Up/(r_L)|^^^, = 

and have units Xj^fA [55]. Two convenient guided-mode 
bases are the quasilinear and quasicircular polarization 
modes, described in Appendix A [54]. 

We consider nanofibers that support only the lowest 
HEii guided modes at the relevant frequency wq [56], and 
thus we drop the mode index j. In this case there are four 
guided modes: two polarizations p, each with propaga¬ 
tion constant P{uio) = ±/3o corresponding to forward and 
backward propagation. The guided-mode contribution to 
the dyadic Greens function is then 


Gg(r, r'; wq) = / d/3 


- 2 wg 

-00 ^wg-a; 2 (;/ 


-U/3,p(rJu^,p(ri)e 


il3{z—z') 


( 8 ) 


1 

(uJo+ief-uj'^iP) 


2w(/3) 


uo+ie-uj{P) 


1 

wo-l-ie-l-w(/3) _ 


and employ the usual distribution identities [42], 


+iTr6{oJoToj{P))- ( 10 ) 

Only the positive-frequency component contributes to 
the ^-function, and it follows that the imaginary part 
of the Green’s function at r = r' that determines the res¬ 
onant Purcell enhancement of spontaneous emission into 
the guided modes is [41, 45, 57] 

Im[G^+\r;r';a;o=Wep)] =7r‘^^ut,p(ri)u^p(r)_), (11) 

c>,p 


lim -^-— = V 

^ 0 + UJod-ie^UJyp) 


where uj^g is resonance frequency of the atomic scatterer. 
The energy level shift of the scatterer due to its proximity 
to the dielectric is found from the real part of the Green’s 
function at r = r'. To find the total modified sponta¬ 
neous emission rate and energy level shift one must in¬ 
clude the unguided radiation modes [49] or employ other 
representations of the Green’s function [43]. 

Equation (9) is the central result from which we can 
calculate the dispersive response. Gonsider a forward- 
propagating input field in the guided modes with fre¬ 
quency Wo, positive-frequency amplitude and ar¬ 

bitrary polarization, E|j|’^(r) = J'g''’^Uin(rj_)e*^°^ disper- 
sively coupled to an atom at position r'. The effective 
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mode area at the atom’s position is determined from 
the total cycle-averaged power transported along the 
nanofiber, = (vg/2'K) J n^(r_L)|E|j|^(r)p, and 

the intensity at the atom, lin(r') = (c/27r)|E|j^^(r')p, 
via the relation [58], 

^ \ (12) 

/i„(r') rig I Uin(r(L)P 

where rig = c/vg is the group index of refraction. 

Substitution of the guided-mode Green’s function, 
Eq. (9), into the Lippman-Schwinger equation, Eq. (2b), 
yields the transmitted (forward-scattered) and re¬ 
flected (backward-scattered) output fields, Eout(r) = 

[u+,out(i’±)e*^“^ -f u__out(r±)e"*^‘’^], 

U-H ,out Cptppf u+,p/(r_L) (13a) 

,out (■■x) = E Cpr’pp'u_ p'(rj_), (13b) 

p,p' 

where we have decomposed the input into the polariza¬ 
tion eigenmodes, Uin(rj^) = Z)p CpU+,p(rj_). For z > z', 
the transmission and reflection matrices are 

tpp' =Sp,p' -f 2TTikong u)]. ^(r^L) • a ■ u+,p/ (r^L), (14a) 

Tpp, =2TTikong u)l_p(r(L) • a • u+,p/ , (14b) 

We focus here on the transmitted fields whose interfer¬ 
ence with the input field for z > z' results in a phase shift 
and a polarization transformation. For weak scattering 
the diagonal terms, tpp ~ -y/l — determine the 

phase shift and attenuation induced on each polarization 
mode, 


5(j)p = ^^Re(app), (15a) 

^in 

= ^^lm(app). (15b) 

^in 

Here, the {p,p'}-element of the tensor polarizability is 
given by, appi = ep ■ a ■ ep, with unit vectors for 
each of the forward-propagating mode functions, ep = 
u-e.p(i'l)/|u+,p(r(L)|. 

The phase shift per atom, Eq. (15a), is modified over 
free space in two ways, both of which are captured 
by the effective mode area Hin. First, although mate¬ 
rial dispersion in an optical fiber is negligible over the 
distances we consider, additional waveguide dispersion 
can lead to a significant reduction in the group veloc¬ 
ity [8, 14]. Such “slow light” enhances the atom-photon 
coupling strength. In the nanofiber geometry this ef¬ 
fect is moderate - we calculated the group index to be 
rig « 1.40. Second and more importantly, the tight spa¬ 
tial confinement as measured by OD /Na significantly in¬ 
creases the coupling strength over free space for every 
atom along the nanofiber, which yields strong coopera- 
tivity. In contrast, in free space diffraction restricts the 


collective phase shift for an ensemble of atoms [15, 59]. 
For a Gaussian beam with beam waist wq, the total phase 
shift induced by a collection of polarizable atoms will be 
St/> = Nes 27 rkoRe(a) /A, where A = 7riyg/2 is the beam 
area at the focus and Neff is the effective number of atoms 
that radiate into this mode. One can couple strongly to 
few atoms at the center by tightly focusing the beam or 
couple weakly to many atoms by choosing a larger focal 
volume, but hence, smaller cooperativity per atom. 

The off-diagonal terms in the transmission matrix, 
Eq. (14a), describe the polarization transformation. For 
example, if we take the polarization of the modes to be 
the quasilinear, p = {H,V} as defined in Eq. (A4), then 
tnv = XFar is the rotation angle of the Stokes vector on 
the Poincare sphere corresponding to the Earaday effect 
[60, 61]. The phase difference in that basis, 5(j)H — S(j)v, 
corresponds to birefringence induced on the guided mode 
and tfiv to Earaday rotation. Analyzed in the quasicir¬ 
cular polarization modes (p = ±), given in Eq. (A3), the 
differential phase S(j)+ — 6(j)- corresponds to Faraday ro¬ 
tation and _to birefringence. We make use of such po¬ 

larization transformations as a means to nondestructively 
measure the atoms and generate collective spin squeez¬ 
ing. 


III. HEISENBERG-LANGEVIN-PICTURE 
SOLUTION AND ATOMIC RESPONSE 


The Lippmann-Schwinger solution, Eq. (2b), deter¬ 
mines the input-output relation for linear atomic re¬ 
sponse given by the polarizability tensor a. In this sec¬ 
tion we connect this with the fully quantum mechani¬ 
cal description of dispersive atomic response and input- 
output relations for the quantized guided modes. Follow¬ 
ing Ref. [49] , we use a Heisenberg-Langevin approach for 
one-dimensional systems. 

The positive frequency component of the quantized 
electric field operator decomposes into guided and ra¬ 
diation (unguided) modes, E^+) = Eg"''^ -f Er~''\ where 



(16a) 

(16b) 


The HEii guided modes are specified by p. = (a;,6,p), 
where uj is the mode frequency, p is the polarization, 
and the propagation direction 6 = ± corresponds to 
wavenumber hf^ipj). The radiation modes are specified 
by = {uj, j3,m,p), where m is the azimuthal (an¬ 
gular momentum) quantum number, p labels the two 
orthogonal polarizations, and longitudinal propagation 
constant /3 can vary continuously from —kn 2 to kn 2 , 
with k = uj/c [42, 49]. The creation/annihilation op¬ 
erators satisfy the usual continuous-mode commutation 
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relations, = 6b.b'Sp^p'6{u} — co') and [a,y,a]^,] = 

Sm,m'Sp,p'S(uJ - oj') 6{(3 - 13 '). 

The Hamiltonian for the system is 


H = HF+HA+Hinu (17) 


where the free-held Hamiltonian decomposes into guided 
and unguided modes, 



Integrating the held equations. 



(24a) 


+ *E5:,e., /‘dt'e—(24b) 

JU 


We consider here alkali atoms with ground and excited 
levels, {\g) = |nS'i/ 2 ,/, to/)}, {|e) = |nP//,/', to//)}, 
where |/, to/) denotes the hyperhne sublevels. The free 
atomic Hamiltonian is 

Ha = '^ Egdgg + ^ Eee , (l9) 

9 e 

where aij = |t)(j|. In the rotating wave approximation, 
the atom-held interaction Hamiltonian is 


substituting into Eq. (23c), and making the usual Markov 
approximation [49] gives an expression for the ground- 
excited coherences. This yields 


ddne .. V—^Egg/ ^ ^^ 

■ = -lUJeg(Tg^-2_^—(Jge' +%2_^ 


dt 


dw (to) +E / 

^ b,p m,p ~ 


{,^e,e' ^gg' ^g,g' ^e'e) 
eig 
kn2 

dl3g^^e',g'd^{to) |.e““ 
'—kn2 


Hint=-dE = - 




[dgg.E(+)(r')-f4e-E(-)(r')] , (20) 


where the atomic dipole operator is projected between 
excited and ground subspaces, dgg = PedPg. The inter¬ 
action Hamiltonian then takes the form. 


Hint = “ E E / ^99,e,g dp deg 

e,g y b,p ® 

poo pkn2 \ 

+E / i^,e,g ^eg + H.C., (21) 

m,p-70 J-kn2 J 

where the coupling constants for guided/radiation modes 
are 


dgp,e,g = \\— (e|d|g) • u^(r(L)e*'’'^(‘^)^ 

V ^3 

dgy,e,g = (e|d|ff) • u^(r(L)e*^(‘^)"^. 
The Heisenberg equations of motion are 

''Y.9l,e,g 


ddp 


dt 

= —icud^ 

ddi, 


dt 

= —icudi 

d&ge 

— ^^65 < 

dt 

{ 

E 99,e', 


b,p 




e,9 


d-E 


pkn2 


' —kn2 


(22a) 

(22b) 

(23a) 

(23b) 


ipee'dgg' Sggide'e) 


(23c) 


where the decay rates of excited-populations and coher¬ 
ences are given by 

Eee' = 27 r ^ ( gp.e.gffp.e'.q 

9,9 

/ kn2 

df3 gi/^e,ggi/^e' ,g\‘^=t^cg^ (25) 

and the small energy shift is absorbed into the transition 
frequency Wgg = {Ee — Eg)/fi. Equation (25) captures the 
modihcation of the spontaneous emission rate due to the 
nanohber. The hrst sum describes decay into the guided 
modes and the second into the unguided radiation modes 
[43, 49, 62-64]. The decay rate of a given excited state 
into all guided modes is given by 

re°=27r^|g;x,e.glL..,=^—El(®|d|5')-Uhp(rD|^. (26) 

L ^ ^g L 

o,p,g ^ b,p,g 


This is in agreement with the expected expression from 
the guided-mode contribution to the dyadic Green’s func¬ 
tion in Eq. (II), 


Eg^ = ‘^'= 9 ) -(eldlg), (27) 


which is enhanced over the free-space rate by the Purcell 
factor. 

Here we are interested in linear response for excitation 
far from resonance. We follow Ref. [65] and consider an 
atom sufficiently far from the fiber surface such that the 
modification of the spontaneous emission rate is small. In 
this case the decay rate is approximated as Egg/ « Se,eHe, 
where Eg is the total decay rate from excited state |e), 
given by the diagonal elements of Eq. (25). In steady 
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state, the dipole operator in the linear regime (&ee' —t 0) 
is approximately 


'ge 




9' 


9^,e,g' 


i'kn2 


, -- --eg 

o,p ^ 


-\-iTe/2 


a/x(io) 


+E/ d/3. _ 

J—kn2 ^ Weg“l“^Ie/2 


.(to) 


m,p 




• (28) 


By substituting this into Eq. (24a) and defining asymp¬ 
totic modes, a“(a;) = limtQ_j._oo a(to)e*“*“, a°“*(a;) = 
limt_>+oo a(t)e*“‘ [51], we obtain the input-output rela¬ 
tionship for the guided modes. 


t / \ ^in/ \ o ■ \ \ ' ^ 9g,,e,g3 ^in / \ 

(w) = (w) -27rz 2^ 2^ ^gg' . , ,t. /o (‘^) 


V ,p' e,g,g' 


^ ^eg 3“ tEg/2 


-2^*E E 

m,p e,g.g'' 


(.fcrio 


/ J ~ 9g,,e,g9v',e,g' , /■r)n\ 

/ d/3ogg.———— a, (w). (29) 

'— kri'i ^ <-^eg~TlLel ^ 


This input-output relation contains the phase shift on 
forward scattered modes as well as attenuation due to 
elastic scattering into all other modes. For a probe with 
frequency ojo, Eq. (29) agrees with the expected form 
given by the Lippmann-Schwinger equation in the first 
Born approximation [66], 

Eiut,g(i’,wo) = E[^g(r,u;o) 

+ (r, r', Wo) -a- [E|+^(r', a;o)-KE^^(r', wq)] , (30) 

by noting that for the guided-mode dyadic Green’s func¬ 
tion given in Eq. (9), 


J u*(rw) • G^+)(r,r',wo) • 



. 27rwo ^ ^ ^ 

=z-Uft p(rj_) • a • Ufc-p/(r^). 


(31) 


Here, the atomic polarizability operator [61, 67, 68], is 
given by 


A 1 \ (5|d|e)(e|d|5') , ,, 

“ = -T E lg) A , .T. /o (g I 




■ /rg/2 


(32) 


and Agg = wq — uieg is the laser detuning from the atomic 
transition. For an atom in ground state \g) and polariza¬ 
tion p, the phase shift can be expressed as [65] 


^^p,g 


Stt —u2p(r(L) • Re[(5|a|5)] • u+,p(r'L) 
^g 

wq 27r|(e|d|g) ■ u+,p(r'|_)p 
Ug ^ hAeg 


We employ this dispersive response for QND measure¬ 
ment of atoms, as we describe in the next section. 



FIG. 2. Quantum interface for spin-polarization coupling of 
two ID lattices of cold, trapped atoms and the guided modes 
of an optical nanofiber, a) Schematic of the interface. A lin¬ 
early polarized probe is launched into the nanofiber and the 
output light is analyzed in a polarimeter. The atoms (green 
circles), trapped in the x-z plane, couple to the evanescent 
portion of the guided H- and E-modes. Contours of the H- 
and E-mode intensities in b) the ^-direction and c) the trans¬ 
verse x-y plane show the mode anisotropy at the atomic po¬ 
sitions. 


IV. QND MEASUREMENT OF ATOMS 

The dispersive interface between the atoms and 
nanofiber guided photons provides the entangling mech¬ 
anism necessary to perform a QND measurement on 
the atoms. We restrict here to the quasilinear modes, 
p = {H,V}, of a single HEn guided mode at frequency 
Wo, whose form is given explicitly in Eq. (A4). In typical 
experimental configurations, two one-dimensional arrays 
of atoms are trapped on either side of the nanofiber, see 
Fig. 2. We define coordinate axes {x,y,z) with z ori¬ 
ented along the fiber axis for forward propagation, and 
the two chains of atoms lie in the x-z plane at azimuthal 
angles </)' = {0, tt}. In the evanescent region, the iJ-mode 
is purely e 2 ,-polarized aX (j) = -Xn:l2 and the E-mode is 
purely Gy-polarized at (j) = {0,7r}. At other azimuthal 
angles the electric field is generally rotating along an el¬ 
lipse in the x-z plane. The atoms at (j)' = t) experience 
H and V fields, 

Ub,H(r±,(p = 0) =V2[e^Ur(r±) -I-ibe^u^(r±)] (34a) 

Ub,v(r±,(/) = 0) =V2eyU^{r±), (34b) 

where the real-valued functions Ua{r±), given in 
Eq. (A5), depend only on the radial coordinate. On the 
opposite side of the fiber at (j)' = atoms experience 
the same transverse electric field, but the z-component 
changes sign. This broken symmetry has been used to 
selectively address and separately control the two atomic 
arrays [30, 32, 36]. 

We consider quasi-monochromatic fields at carrier fre¬ 
quency Wo that are sufficiently narrowband. Aw ^ wq. 
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For each guided mode we define input propagating, 
continuous-mode field operators in the interaction pic¬ 
ture [46, 47, 50], 


1 

ab,p{z,t) = 

V Stt Jo 


(35) 


ear HV-hasis, 

Si{t) = 5 [a\j{t)aHit) - alr{t)av{t )], 
S2{t) = 5 [alj{t)av{t) + al{t)dH{t)] , 
-^3 W = ii [«ff (Oay(i) - a|,(t)aff (t)], 


(40a) 

(40b) 

(40c) 


that satisfy the free-field commutation relations, 

[db,p{z,t),dl,p,{z',t')] = Sb,b'Sp,p'6{t-t-{z-z')/vg). (36) 

In terms of these propagating modes the quantized elec¬ 
tric field operator, Eq. (16a), becomes 

(rj_, (j), z; t) =y^W ^^^^° Ub,p(rj_, 4>)hA^ (37) 

h,p V 

Considering here only the forward-propagating guided 
modes (6 = -I-), we drop the b index. The propagating 
electric field, Eq. (37), interacts with the trapped atoms 
via the dispersive light-shift Hamiltonian [61, 68, 69], 

Na 

Hp^s = -Y. (4; t), (38) 

n—1 


where is the atomic tensor polarizability operator, 
given in Eq. (32), for the atom trapped near the 
nanofiber surface at position r'^. We ignore here any 
effects of atomic motion and treat the atoms as localized 
at fixed positions in space. 

The Lippmann-Schwinger scattering equation, 
Eq. (30), follows in the time domain as the evolution 
of coarse-grained input-ouput modes [46, 51, 65]. Since 
multiple scattering is negligible and the propagation 
time across the ensemble is small compared to the 
atomic dynamics, we drop the position label z and index 
the propagating fields by time alone; db^p{z,t) —>■ db,p(t) 
as is standard in input-output theory [46, 70] . It 
follows that the effects of retardation can be ignored, in 
which case each term in the sum over atoms contributes 
equally for all atoms (for details see [50, 69]). The 
forward-propagating output fields are then given by the 
Fourier transform of Eq. (29), yielding [50] 

a°"*(t) = a”(t) [A^oU*(r(L,0) • a • Up. (r^L, 0) 

zj p/ 

+ N^u*p{r'j_,TT) ■ a ■ Up>{r'j_,Tr) 


a”(t), (39) 


where {NqjNtt} are the total number of atoms trapped 
at (/)' = {OjTt}. The quantum effects from the first term 
give rise to shot noise in the transmitted field at the de¬ 
tector. The second term represents scattering into the 
guided modes, as described by the dyadic Green’s func¬ 
tion, Eq. (27). 

We introduce the vector Stokes operators that describe 
the polarization of the propagating fields in the quasilin- 


that satisfy equal-z commutation relations following from 
Eq. (36), 


=ieijkS{t-t')Sk{t). (41) 

These, along with the total photon flux operator, 

So{t) = i [dlp{t)dHit) + alr{t)dv{t )], (42) 

are used to reexpress the Hamiltonian, Eq. (38), in the 
HV-hasis, 


Hls = - ^TThkoUg ^ N^> 

tp' — 0,7T 

I [KHHi4>') + Kvv{(l>')] So{t) 

+ [Khh{(I>') - Kvv(,4'')] Si{t) 

+ ((/)') -b iFvff (<(>')] <^ 2(0 

+ i[kHv{^')-KvH{cl^')]S3{t)y (43) 

The atomic couplings to the {H,V} modes, 

kpp'ik) = |u;(r(L,</>')||up'(r-(L,</>')! app/((/)'), (44) 

are determined by components of the quantum mechan¬ 
ical tensor operator weighted by the transverse mode 
functions at the atomic position, appi{(j)') = e*,{(j)') ■ a ■ 
ep((j)'), whose classical analog appeared in Eq. (15a). 

We explore a QND measurement of ^^^Cs atoms in 
the electronic ground state, 65'i/2, via polarization spec¬ 
troscopy based on the collective atom-light coupling 
described by the dispersive light-shift Hamiltonian in 
Eq. (43). Polarization transformations occur due to the 
tensor nature of the atomic response, 

^ = (45) 

/,/' ij 


v-r /s 

where the operator A(/,/') = kjif, f')e^ 0 
decomposes into irreducible components within each 
ground hyperfine multiplet / for light detuned near ex¬ 
cited multiplet /', 




(1). 


+c' 




(46) 


Here, ao(A//-) = -^ A„,+ir /2 is the characteristic 
dynamic polarizability where tJo = 3A^/27r is the resonant 
scattering cross section, f is the atomic spin operator 










in hyperfine multiplet /, and C^/f) are coefficients for 
irreducible rank-X components defined in [61]. 

In addition to the atomic tensor response, the 
nanofiber geometry gives rise to unique features of po¬ 
larization spectroscopy not present in free space. The 
spatial anisotropy of the intensity for the quasilinearly 
polarized guided modes leads to unequal scattering of 
the H and V modes, producing intrinsic birefringence 
even for a purely scalar atomic polarizability. In partic¬ 
ular, atoms trapped on the quasi-i/ axis leads to a phase 
delay of this mode relative to the fast quasi-y axis. This 
birefringence was exploited by Dawkins et al. [.33] as 
a mechanism for implementing a dispersive QND mea¬ 
surement of the number of atoms trapped around the 
nanofiber, as we treat in the next section. 

A. Dispersive atom number measurement 


integrated measurement is described by the operator 
jCi = The shot-noise variance of the po- 

larimeter, AAI^jsN = Xn^lT for integration time T, 
determines the fundamental resolution of the polarime- 
ter. The smallest detectable atom number using this dis¬ 
persive measurement is thus, 5Na [71]. 

In an ideal setting, 5Na can always be reduced by in¬ 
creasing the integration time, but in practice this time is 
limited by atom loss. As a coarse approximation we take 
this time to be T = 77^, where 7 s is the photon scatter¬ 
ing rate in free space, and assume perfect quantum effi¬ 
ciency of the detectors. For detuning A large compared 
to the excited hyperfine splitting on the Dl- or D2-line 
{f = 1/2 or 3/2), the unit-oscillator scattering rate is 
7 s = ^ with effective area determined by 

the probe at the atomic position, Eq. (12). In this limit, 

the rotation angle xn = Cj°^(o’o/^Ar)(r/2A), Eq. (48), 
yields a shot noise-limited atom number resolution. 


The anisotropy of the guided modes provides a mech¬ 
anism for counting the number of atoms trapped around 
the nanofiber based on polarization spectroscopy. We 
consider Na atoms, each in a completely mixed hyperfine 
spin state. In this case the atomic polarizability tensor 
in Eq. (46) reduces to (Ay(/,/')) = and the 

collective interaction is determined entirely by the the 
scalar (rank-0) terms. With the atoms trapped along the 
quasi-axis, while (Khh) ^ {Kvv), the off-diagonal 
elements in Eq. (43) do not contribute to the Birefrin- 
gent interaction we are interested in and actually vanish 
{{Khv) = (Kvh) = 0) when x—, y— or z—axis is chosen 
as the quantization axis which includes the one close to 
the optimal choice of quantization axis for spin squeezing 
we will discuss in the next section. Atoms on either side 
of the nanofiber experience the same scalar light shift 
yielding from Eq. (43) the Hamiltonian for QND mea¬ 
surement of atom number, 

Hn = —2'Khkong E (</')) - {kvvi^'))]Si{t) 

= hxNNAk(t). (47) 


This birefringent interaction induces a rotation of the 
Stokes vector around the S'l-axis on the Poincare sphere 
through an angle. 


Xn 


^0 r^(o) r 

^N “ 2Aff, ’ 


(48) 


characterized by an effective area, = 

/2) (IUff (ri)l2 -1uy (ri) 12). 

Dawkins et al. [33] used this interaction to make 
a dispersive measurement of Na via birefringence po- 
larimetry in the usual way: launching linearly polarized 
light at 45° to the quasi-axis, Uin = {uh + uy)/-\/2, 
and measuring the differential power between the guided 
right-and left-circularly polarized photons. Thus, the 


6Na 


(7j/) y AinCro ’ 


(49) 


where = J2f f the far-detuned, rank-0 co¬ 

efficients on a j — j' transition [61]. Using the param¬ 
eters reported by Dawkins et al. [33], we find the shot- 
noise limited minimum detectable atom-number 6Na 
10 for atoms trapped at 1.8a ^ 2.0a from the fiber axis 
with a D2-line probe light. 

In practice, loss and decoherence limit the atom- 
number resolution [33, 72]. The experimental implemen¬ 
tation reported by Dawkins et al. [33] implies a resolution 
of a few tens of atoms for 200 ~ 1000 trapped atoms. A 
similar experiment based on a two-color QND measure¬ 
ment in a nanofiber geometry was recently carried out 
by Beguin et al. [34] to squeeze the uncertainty in the 
number of trapped atoms. They achieved an atom num¬ 
ber uncertainty of 6Na = 8 for Na ~ 2500 atoms, well 
below standard quantum limit, 6Na = ^/Na- 


B. Collective spin squeezing via QND 
measurement 


The same birefringent interaction, Eq. (43), can be 
utilized in a QND measurement to squeeze the pro¬ 
jection noise of the collective atomic spin. We con¬ 
sider squeezing of the uncertainty associated with the 
“clock states” of cesium, jf) = l65'i/2, / = 4, m/= 0) 
and jl) = 16 iS'i/ 2,/ = 3,m/ = 0), which define a pseu¬ 
dospin within each atom and associated Pauli operators 
{(Ti, CT 2 , da}. The quantum uncertainty in the collective 
pseudospin, 

, Na 

= (50) 

n—1 
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fundamentally limits the precision of atomic clocks [73]. 
For atoms prepared in a spin coherent state (SCS) the 
projection noise, AJ^\^^^=Na/4:, sets the standard quan¬ 
tum limit for spin measurements. A spin squeezed state 
(SSS) exhibits reduced fluctuations, AJ||ggg< A 4 / 4 , due 
to negative pairwise correlations between the atoms [74] . 
Spin squeezing is typically quantified with the metrolog¬ 
ical squeezing parameter defined by Wineland et al. [73] , 

A P 

^2 = ^ (51) 

(J||)2’ 

where (J||) is the mean collective spin along the direction 
of spin polarization. 

The clock states are defined to have zero projection of 
angular momentum with respect to a bias magnetic held 
that dehnes a quantization axis, e^. Within the clock- 
state subspace the rank-1 vector light shift in the dis¬ 
persive Hamiltonian, Eq. (43), vanishes since (t|/fc|t) = 
('I'l/fcj'l') = 0 for any direction of the spin, k, and any 
quantization axis, ej. Furthermore, as shown below, 
atoms on either side of the nanohber experience the same 
birefringent coupling. The resulting Hamiltonian, re¬ 
stricted to the clock subspace, couples the guided held 
of the nanohber to the Ja-component of the collective 
pseudospin. The interaction has contributions from both 
the scalar and tensor light shifts, 

+ XV,t) - ixH,i + XV,i)] JsSoit) (52) 
+ [{xH,-r - xvt) - {xH,i - XEt)]'^3<Si(t)|, 

where the coupling strength between an atom in the clock 
subspace and a photon with polarization p = {H,V} is 

Xp,f = -2TTkong\up{rj_)\‘^{f,0\app\f,0), (53) 

and / = {4,3} labels The diagonal terms in the 

polarizability tensor are the same for atoms at positions 
above and below the nanohber, and thus all atoms con¬ 
tribute equally. In addition, a constant birefringence pro¬ 
portional to JqSi is neglected here as it can be canceled 
with a compensating waveplate. Finally, the hrst term 
in Eq. (52) does not affect polarization spectroscopy, but 
will act to rotate the pseudo-spin around the Ja-axis 
of the generalized Bloch sphere proportional to classi¬ 
cal intensity huctuations. While this does not affect the 
squeezing of projection noise in J 3 , it affects the metro- 
logically relevant squeezing by adding uncertainty to the 
direction of the mean spin. By choosing a “magic fre¬ 
quency” at which the light shifts on the two clock states 
are equal, 

XH,t + Xy.T = Xh,i + Xfxi (54) 

this term can be canceled [75], where we have ignored the 
imaginary part of the coupling strengths in the dispersive 
regime. Using the Dl-line of ^^^Cs atoms as the probe 



FIG. 3. Parameters of the atom-light interface using the 
clock states of ^®®Cs. (a) Energy level structure for atoms 
probed with one of two magic frequencies on the Dl-line. (b) 
Magnitude of the magic detunings at which the clock states 
are equally light-shifted, |A/|/27r = jtJo ~ (in units 

of MHz). There are two solutions shown as blue (dashed) 
and red lines, (c) The coupling strength, as measured by the 
magnitude of the polarization rotation angle on the Poincare 
sphere, jXjji) for an atom trapped in the x-z plane at a dis¬ 
tance r'±_ = 1.8a from the hber center. In both (b) and (c) we 
plot the parameter as the direction of the clock-state quan¬ 
tization axis is varied in the x-y plane, for the two possible 
choices of magic detunings. See text for details. 

light, there are two magic-frequency solutions, W3 and 0)4, 
shown in Fig. 3(a). 

Because the guided probe light at the position of the 
atom will generally be elliptical, the light-shift interac¬ 
tion coherently couples different magnetic sublevels in a 
given manifold /, and thus does not conserve J3. For 
example, the ellipticity of the probe light leads to a hcti- 
tious magnetic held proportional to zE|~^(r') x E|j^^(r') 
that causes a precession of the spin within hyperfine man¬ 
ifold /. This can be mitigated by a sufficiently strong bias 
magnetic field compared to the hctitious held [76]. 

The remaining QND interaction Hamiltonian is 

Hja = hxj^JsSiit), (55) 

where the rotation angle on the Poincare sphere at the 
magic wavelength is, 

Xj3 = ixH,t-Xv,t)-ixH,i-Xv,i) = 2{xH,t-XH,i)- (56) 

In the standard way, squeezing the uncertainty in J 3 by 
QND measurement can be generated by preparing the 
atoms in a SCS along Ji, passing a probe prepared along 
S 2 with photon hux TV/,, and continuously monitoring the 
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S'a-component of the guided light in a polarimeter. The 
measurement strength, 

K=\xjfNL, (57) 

quantifies the rate at which we squeeze projection noise, 
with Xjs given in Eq. (56). In the absence of any deco¬ 
herence, such a QND measurement for integration time 
T squeezes the initial uncertainty in J 3 according to 
(AJ|) out — (AJ|)in/(l -hr), where 

r = kT(A J|)i„ (58) 

is the integrated measurement strength [15, 60]. 

The strength of the birefringent interaction arises from 
two fundamental sources. The anisotropy of the H and 
V polarized modes leads to a polarization-dependent in¬ 
dex of refraction, as described in Sec. IV A. In addition 
there is a dependence of the atom-photon coupling on the 
internal spin state of the atom due to the atomic tensor 
polarizability. In particular, we are interested in the de¬ 
pendence on the two clock states of the atom. This spin- 
dependent coupling will depend on the choice of quan¬ 
tization axis that defines the clock state with projection 
rrif = 0 . 

We combine these two effects and obtain a compact 
expression for the coupling strength Xj 3 using the irre¬ 
ducible tensor decomposition of the atomic polarizability, 
Eq. (46). Let {ex,ey,ez} be a space-fixed Cartesian co¬ 
ordinate system, where defines the quantization axis 
of the atom, set by the magnetic held. Because of the 
azimuthal symmetry of clock state around the e^-axis, 
the polarizability tensor is diagonal in that basis. Noting 
that (/,0|A|/,0) = 0 and (/,0|/||/,0) = (/,0|/||/,0) = 

(/, Ojf^l/, 0)/2 = /(/ -I- l)/2, it follows that the expec¬ 
tation value of the irreducible rank -2 component of the 
atomic polarizability is 


/' 

The combined scalar and tensor light shifts yield a cou¬ 
pling strength, Eq. (53), 

XpJ = ngCro{af \up{r'j_)f - bf |e^ • Up(r(L)|^), (60) 

with coefficients that depend on detunings and atomic 
structure. 


which depends on the choice of quantization axis. 

We write the effective rotation angle in the Hamilto¬ 
nian, Eq. (55), as 


with an “effective detuning” set by the magic-wavelength 
condition. 


^73^ - f 




and an effective area given by 

1 _ I • u V (ri)| 2 (ri)| ^ -1 • u// (ri)| ^ (ri)| ^ 

^ J 3 “ ”9 


|Uff(ri )|2 -H |uy(ri )|2 


( 66 ) 


We see here the explicit dependence of the coupling 
strength on both the anisotropy of the modes and on the 
tensor atomic response, which in turn depends on a par¬ 
ticular choice of clock states. The quantization axis that 
maximizes Xjs is that which minimizes Aj^ at a given 
magic detuning. Since the ^-component of the guided 
modes is 90° out-of-phase with the transverse compo¬ 
nents, the quantization axis maximizing the atom-light 
coupling is specified by an angle in the transverse x-y 
plane, (p, 


ej = cos (pBx + sin pey. (67) 

The dependence of the magic detunings on the direc¬ 
tion of quantization axis is shown in Eig. 3(b) for atoms 
trapped at a typical distance of r'j_ = 1 . 8 a on the x- 
axis. In typical operating regimes, the magic frequencies 
are hundreds of MHz from resonance with either excited 
state, placing the interaction in the off-resonant, disper¬ 
sive regime. Using these magic detunings, in Fig. 3(c) we 
show the variation in Xj^ a function of (p. This sug¬ 
gests that, based solely on the strength of the coherent 
interaction, the cc-axis is the optimal quantization axis. 
As we will see in the next section, the optimal quantiza¬ 
tion axis is significantly modified when decoherence due 
to optimal pumping is included. 


C. Decoherence due to optical pumping 


p-iio) 


“/ = V ( 

J o 


/(/ + ^) ^( 2 ) \ r 

6 ^) ‘iAffi ’ 


/' 


4A 


//' 


At the magic wavelength set by Eq. (54), 

aj - «3 _ \ei ■ uy(r(|_)|2 -q [eg ■ ug(r(i_)|2 

54-53 |uff(r'i _)|2 -t |uy(r'j _)|2 


(61) 

(62) 


(63) 


The treatment above considers an idealized QND in¬ 
teraction. The coupling of the atoms to the probe, how¬ 
ever, will always lead to scattering of photons into modes 
other than the forward-scattered guided mode. This is 
accompanied by optical pumping that destroys the en¬ 
tanglement associated with spin squeezing. In addition 
it reduces the metrologically useful signal. The maxi¬ 
mum achievable metrologically relevant squeezing is de¬ 
termined by the balance of this decoherence with the 
QND measurement. 
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We model this using a first-principles stochastic master 
equation description (SME) [15, 77], 

dp = sJ'^T-L[p]dW + '^C[p]dt -t ^I?„[/5]dt, ( 68 ) 

' n 

where s = sign(xj 3 ) and p is the collective atomic state. 
The measurement strength n = determines the 

rate of the spin squeezing in the absence of decoher¬ 
ence. The first two terms describe the QND measure¬ 
ment, where dW is a stochastic Weiner increment satis¬ 
fying dW^ = dt. The conditional dynamics that result 
from the measurement are generated by the superopera¬ 
tor 


We calculate the squeezing parameter as a function of 
time based on the evolution of atomic correlation func¬ 
tions, where operators evolve according to the adjoint 
form of the SME in Eq. ( 68 ) . The collective atomic vari¬ 
ables obey the following stochastic equations of motion 
(see Appendix C), 

dNc = -jooNcdt + 27o3(J3)dt (73a) 

d(Ji) = -7ii(Ji)dt (73b) 

d(J3) = sy/nAJ^dW - 733(^3)^^ + ^jsoNcdt (73c) 

dAj| = —«;(Aj|)^df — 2733A j|dt 

+ 2 (2733-700) (703-2730) (i 3 )di, ( 73 d) 


H[p]=J3P + pJ3-2{J3)p, (69) 

and the collective Lindblad map is 

£[/!]=- i {pJi + JIp) + hpJ^. (70) 

The final term in Eq. ( 68 ) describes the effect of optical 
pumping acting locally on each atom along the nanofiber. 

The optical pumping map is governed by a standard 
master equation [61]. Restricting to the two-dimensional 
subspace associated with the clock states, the action on 
the atom is 

'Dn[p\ = E{-y 

/=3.4 

/=3.4 


Here, 7 / is the total rate of photon scattering by atoms 
in state ]/, 0 ) and is the rate of optical pumping be¬ 
tween the clock states, ]/, 0) —)■ l/,0) (see Appendix B). 
Expressed in terms of Pauli operators on the clock-state 
pseudospin, the map acts as 


Tdn[p\ — 


2(7t + 7l) - 7T^T - 71^1 


P 


7t->t + 71-s-l - (n) ^~(n) 
"3 P"3 


, {fi) ^ ^ yiL I , ^ yiL) ^ ^ yn 

+ 7t^icr7 VctY + 7i^to-+ P^- ■ 


(72) 


There are three important features of this map that are 
not typical in a QND measurement of ideal spin-i par¬ 
ticles. Eirst, the map is not trace preserving because 
atoms can be pumped out of the clock states. Second, 
unequal rates of optical pumping for jf) and jj,) polarize 
the mean (J 3 ) towards a value different from that found 
in the QND measurement. Third, owing to the large 
ground hyperfine splitting, photons arising from optical 
pumping of /—>■/ = 3 and / —?► / = 4 are distinguish¬ 
able, thus these processes destroy coherences between jf) 
and ll). 


where the decay and feeding rates are given in Eq. (C4). 
The total number of atoms in the clock-state subspace 
is given by Nc, which primarily decays at rate 700 . The 
final term in Eq. (73d), proportional to (J 3 ), is typically 
negligible since in most applications (J 3 ) ^ Nc- We re¬ 
tain this small correction since unbalanced optical pump¬ 
ing acts to polarize the atoms and alters the rate of atom 
loss. 

To hnd the peak squeezing in the presence of optical 
pumping, we numerically integrate Eqs. (73a-73d) and 
then use Eq. (51) to calculate the metrological squeezing 
parameter, as a function of time. We choose here the 
magic frequency close to the / = 4 /' = 4 transition, 

d; 4 , which is furthest from resonance with both excited 
hyperfine transitions. Typical time evolution is shown in 
Fig. 4 for 2500 atoms trapped a distance r'± = 1.8a from 
the center of the nanohber, where time is scaled to the 
characteristic scattering rate. 


Is = 


ao r2 
7lin4A2^ 


Nl- 


(74) 


We study the dynamics for two choices of quantization 
axis: (i) along the x-axis and (ii) along the numerically 
determined optimal axis. Figure 4(a) shows the time evo¬ 
lution of the squeezing parameter. We achieve a maxi¬ 
mum squeezing of 4.7 dB when the clock states are chosen 
along the optimal axis; ipopt « 86 ° in Eq. (67). 

The peak squeezing is ultimately limited by the com¬ 
bined effects of optical pumping on both (Ji) and AJ|. 
Here, as in a free-space model [15], the primary factor 
that limits metrological squeezing is the decay of the col¬ 
lective mean spin (Ji). A scattered photon eliminates 
the initial coherence between jf) and jj,) within a single 
atom, thus depolarizing (Ji). Atoms optically pumped 
to magnetic sublevels outside of the clock subspace de¬ 
cay Nc, further reducing (Ji). These effects are captured 
by the depolarization rate 711 in the equation for (Ji), 
Eq. (73b), whose solution is plotted in Fig. 4(b). 

We can gain deeper understanding in the microscopic 
effects of optical pumping on spin squeezing by looking at 
the evolution of the one and two-body correlation func¬ 
tions. In terms of its constituent pseudospins, the collec- 
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FIG. 4. Squeezing on the clock states as a function of time in units of the scattering rate js for 2500 atoms trapped in the 
x-z plane at a distance r'± = 1.8a from the axis of the nanofiber. The optimal quantization axis (red) is compared to the 
quantization along the x-axis (bine). Dashed lines indicate simulations without optical pumping; i.e. no decoherence, (a) 
Metrological spin sqneezing parameter Eq. (51), in dB. (b) Collective mean spin (Ji) and decaying atom number in the 
clock states, Nc (inset), (c) Conditional squeezed variance Ajf. The inset shows the decomposition at the optimal quantization 
axis of Aj| (red dashed) into the single-body variance, fVi(Aj®)^ (green) and the two-body covariance, 

(black), as given by Eq. (75). 


tive variance takes the form 

AJ| = NA(Ajfy + Na{Na - l)(Ajf ) (75) 

for permutationally symmetric states considered here, 
where (1) and (2) label any two atoms in the ensem¬ 
ble. Loss of atoms affects the first (single-body) variance 
term, which scales as Na- The two-body correlations 
which contribute as N'^ to the collective fluctuations, 

(Ajf \ O - (4^V). (76) 

have a much larger influence on the total variance. 
Spin-spin correlations at the heart of spin squeezing, 
rapidly generated by the measurement 
backaction decohere by optical pumping according to 
Eqs. (C2-C3), 


d 

dt 


A.®,. 


) 


op 


- 2733(d®® df'^)+73 o( iq O df+ d%i^^), 


(77) 


where = (|t)(tl + li)(4'l)^”^ is the single-body projec¬ 
tor onto the clock states. In addition, atoms that return 
to the clock subspace after scattering a photon inject ad¬ 
ditional noise into AJ|. All of these effects are included 
in the equation for A J|, Eq. (73d), whose overall and de¬ 
composed dynamical evolutions are shown in Fig. 4(c). 

With our model, we explore optimal conditions for gen¬ 
erating spin squeezing. The choice of quantization axis 
that defines clock states affects both the measurement 
strength and the relative rates of optical pumping. We 
plot the peak squeezing as a function of the direction of 
in the x-y plane in Fig. 5(b). We gain insight into the 
tradeoffs between QND entangling interaction and deco¬ 
herence by independent inspection of the measurement 
strength and optical pumping rates. First, the rate of 
squeezing is determined by the effective optical density 
per atom on resonance, 

OH/Na ^ - = (78) 



FIG. 5. Dependence of the parameters that determine metro¬ 
logical squeezing on the direction of the quantization axis that 
dehnes the clock states, ej. In all cases we consider 2500 
atoms trapped in the x-z plane at distance r'x = 1 . 8 a from 
the axis of the nanofiber. In (b-d) is confined to the x-y 
plane, where the optimal peak squeezing occnrs. (a,b) Peak 
achievable squeezing at the maximum time, measured in dB, 
as a fnnction of the direction of ej. (c) OY>/NA=aoAin/ 

Eq. (78). (d) Rates of atom loss, 700 , and depolarization, 711 
relative to the characteristic scattering rate 7 s. See Eqs. (C4a) 
and (C4e). 


which peaks when ej is along the j/-axis, as seen in 
Fig. 5(c). Choosing along y, the OD/Na is about 
50% larger than along x-axis. The various forms of deco¬ 
herence similarly vary with quantization axis, as seen in 
Fig. 5(d), where we plot the dominant rate of atom loss. 
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7 ooi and the depolarization rate of the mean pseudospin 
(Ji), 7 ii. Because the magic frequency 0)4 is nearly 
equidistant from f = i and /' = 4 when the quanti¬ 
zation axis is near the y-axis (see Fig. 3(b)) this choice of 
quantization axis provides more protection from decoher¬ 
ence. While the decoherence rates in Fig. 5(d) are largest 
near the y-axis, the increase in n more than compensates 
to provide optimal peak squeezing. 

Finally, we explore the optimal conditions as a func¬ 
tion of the trapping geometry. The dispersive entangling 
interaction is based on the collective atomic coupling to 
the evanescent guided-mode fields, which decay exponen¬ 
tially away from the nanofiber surface, as seen in OD/TV^ 
plotted in Fig. 6 (a). From Eq. (73d), the optimal choice 
of quantization axis depends not only on distance from 
the fiber but also weakly on the atom number. Fig. 6 (b) 
because of the competition between squeezing and deco¬ 
herence. At the optimal quantization axis, the strong de¬ 
pendence of peak achievable squeezing on distance from 
the fiber is as seen in Fig. 6 (c) along with the expected 
increase as more atoms contribute to the atom-light in¬ 
terface. 

Several effects limit the reliability of the simulations 
for atoms trapped very near the fiber surface as r±^a. 
First, strong van der Waals interactions modify the light 
shifts and magic frequencies [22, 78]. Second, the optical 
pumping model used here breaks down when the local 
density of states is significantly modified by the presence 
of the dielectric nanofiber [49, 79]. At distances r± > 1.5a 
the atoms’ local environment is roughly that of unmodi¬ 
fied vacuum [49] and a free-space optical pumping model 
in Eq. (72) suffices. 


V. SUMMARY AND OUTLOOK 

We studied the strong cooperativity in the atom-light 
interface that can be achieved based on atoms trapped 
in the evanescent field surrounding an optical nanohber, 
and interacting with a guided mode in the dispersive 
regime. The key parameter that determines the cou¬ 
pling is the resonant optical density per atom. Due to 
the tight confinement of the guided mode over the en¬ 
tire chain of atoms this parameter is OD/iV^ ~ 10“^ 
for typical geometries used in current experiments, which 
approaches that achieved for atomic ensembles trapped 
inside optical cavities of moderate finesse [72, 80]. In con¬ 
trast, the atom-light coupling for atoms in free space is 
typically orders of magnitude smaller, OD/Na ~ 10 “®. 
Under ideal conditions the atom-light interaction is en¬ 
tirely symmetric along the nanohber, providing a plat¬ 
form for long-range correlations independent of distance 
between the atoms. As the light is entirely guided, hber- 
or waveguide-coupled atomic ensembles can be networked 
together or coupled to other physical systems in a hybrid 
platform [81-84] for truly long-range entanglement gen¬ 
eration and distribntion. 

We calculated the dispersive response based on a 


modal decomposition of the dyadic Green’s function, 
which provides a general method to calculate the induced 
phase shifts and polarization rotations of the guided 
modes. With this we studied the QND measurement 
of atoms via polarization spectroscopy. In particular, 
we studied squeezing of the collective pseudospin asso¬ 
ciated the atomic clock states of cesium. The atoms in¬ 
duce a birefringent index of refraction on the light, condi¬ 
tional on the spin state, which provides a mechanism for 
measuring the atomic spin projection and thus squeezing 
its uncertainty. Based on our formalism we calculated 
the nanofiber-enhanced measurement strength that de¬ 
termines the rate of squeezing. 

The peak squeezing one can generate depends on a 
detailed balance between the reduction of spin projec¬ 
tion noise based on QND measurement and the damage 
done to the spin ensemble due to optical pumping. Both 
measurement and optical pumping arise from the same 
physical mechanism - scattering of photons by atoms. 
The former corresponds to cooperative forward scatter¬ 
ing into the guided mode whereas the latter corresponds 
to local scattering into all other modes, primarily the un¬ 
guided “radiation” modes. The cooperativity, specified 
by the effective OD/Na, determines the ratio of these 
two effects and thus the ultimate power of the quantum 
atom-light interface. 

We studied QND measurement using a first-principles 
stochastic master equation model, which allowed us to 
track the atomic correlation functions that define the 
metrologically relevant squeezing parameter. These in¬ 
clude the atomic projection noise uncertainty as well as 
the length of the collective spin vector that defines the 
metrological signal. We find that decoherence acts pri¬ 
marily to depolarize the mean pseudospin and optically 
pump atoms out of the clock subspace, which we treat 
as loss. In addition, optical pumping decoheres the spin 
correlations at the heart of spin squeezing, but at a re¬ 
duced rate compared with the effect on the mean pseu¬ 
dospin. The combined effect of QND measurement and 
decoherence yields a peak squeezing approaching 5 dB 
with ~ 2500 atoms. Larger enhancements in atom-light 
coupling and QND squeezing are possible with modest in¬ 
creases in the number of trapped atoms and/or for atoms 
trapped closer to the nanofiber surface. 

Whereas we have assumed here that atoms can be pre¬ 
pared in a desired clock state defined by a particular 
quantization axis, in practice such preparation will re¬ 
quire optical pumping that may be challenging for atoms 
near the surface of the nanohber. In addition, though 
we have treated the atoms as localized at well dehned 
points, in practice the atoms’ thermal motion can re¬ 
duce the strong coupling described here. Our formalism 
provides a starting point for developing models neces¬ 
sary to study the dynamics of optical pumping, including 
the possibility of cooling atoms to the vibrational ground 
state, where thermal motion is negligible. 

Finally, though we have treated here the case of strong 
coupling due solely to tight conhnement of the guided 


14 




FIG. 6. Parameters that define squeezing as a function of trapping distance and initial atom number Na- (a) OD/Na- 
Inset is the corresponding magic detuning in units of MHz (see text), (b) Optimal quantization axis orientation angle ej in the 
x-y plane for different atom numbers. The line with 500 atoms terminates when the squeezing effect is too weak to be observed 
(r'j_ > 2.3a). (c) Peak metrological squeezing at the optimal ej for different atom numbers. 


mode for atoms near the surface of the nanofiber, we 
can achieve even greater enhancement by combining this 
effect with longitudinal confinement provided by fiber- 
based optical cavities [37, 38, 85-87]. The coupling can 
be further improved under EIT conditions that substan¬ 
tially slow the group velocity [35, 36, 88, 89]. In addi¬ 
tion, quantum control of the internal hyperfine state [90] 
can greatly enhance the entangling power of the atom- 
light interface [91, 92]. For large enough coupling, QND 


measurement should allow production of highly entan¬ 
gled spin states beyond the Gaussian regime [70, 93]. 
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Appendix A: Guided-mode functions for the optical 
nanofiber 


In this Appendix we provide, for reference, the funda¬ 
mental HEii solutions to the homogeneous wave equa¬ 
tion, Eq. (1) with a = 0, for a cylindrical nanofiber of 
radius a and index of refraction given by Eq. (7). At a 
given frequency, uq = ckg, the magnitudes of the longi¬ 
tudinal and transverse wave vectors for a guided mode 
are related by = /3g -b The positive propagation 
constant, /3o = Piuig), is determined from the eigenvalue 
equation that results from enforcing physical boundary 
conditions at the fiber surface [56], 


Jo{ha) 

haJi(ha) 



nf -b 77.2 K'{qa) 1 
2n\ qaKi{qa) h'^a^ 

K\qa) y PI ( I 1 \ 
qaKi{qa)) n\k'^\q^a? K^a^J 


2n 



Inside the nanofiber the transverse wavevector is real, 
= g, where q = y/P q — n^k^, and outside the 
nanohber it is purely imaginary, fc_L = ih, where h = 
\/n\kQ — Pq. The vector eigenfunctions are expressed as 
f;i(r) = (27r)“^/^Uf,_p(rj_)e*^^°^, where the modes are in¬ 


dexed by frequency wq, propagation direction 6 = ±, and 
polarization p. 

A relatively simple form for the guided-mode func¬ 
tions can be expressed in a cylindrical basis {r±,(j),z) 
with longitudinal unit vector , oriented along the fiber 
axis. The transverse unit vectors are related to their fixed 
Cartesian counterparts via the relations 

Gri = Gx cos (/) + By sin (j), (A2a) 

= —Bx sin (j) + By cos (j). (A2b) 

The transverse profile for the quasicircular guided modes, 
p = ±,is 

Ub,±{r±)=[Brjirir±)±ie^u^{rj_)+ibB^u^{rj_)]e^'‘^, (A3) 
and for the quasilinear guided modes, p = {H,V}, is 

Ub.H{r±) =V2[Br^Ur^{r±)cOS(j) 

— B^u^{r±) sm(j)+ibBzUz{r±) cos 0] (A4a) 

U 6 .y(rj.) ='/ 2 [er^Ur^(r_L)sin(/) 

+Btj,u^{r±) cos (/)+* 6 e^U 2 (r_L) sin(()]. (A4b) 

The modes are expressed in terms of real-valued functions 
that depend only on the radial coordinate r±, 

Urj_ir±) =uo[(I - s)Ko{qrj_) -b (I + s)K2{qr±)] (A5a) 

Utj,ir±) =uo [(I - s)Ko{qrj_) - (I -b s)K2{qr±)] (A5b) 

Uz{rj_) (A5c) 

do Ji(na) 

where uq is set by the normalization condition, 
/ d'^r±n{r±)\Ufj,{r±)\‘^ = 1, and Kn are the Bessel 
functions of the first and second kind, f'{x) indicates a 
derivative with respect to the argument x, and 

[J[{ha)/haJi{ha) -b K[{qa)/qaKi{qa)] 

Of particular interest is the z-component, Eq. (A5c), 
which can become appreciable. Note that the phase 
convention in Eqs. (A3-A5) has been chosen to em¬ 
phasize properties of the quasilinear modes and differs 
from that of Le Kien et al. - for instance in Ref. [65]. 
Further details about the guided-mode fields inside the 
nanohber (r_L < a), the radiation (unguided) modes, 
and the quantized form of both can be found in Refs. 
[42, 49, 54, 78, 95]. 


Appendix B: Photon scattering and optical pumping 

rates 

In this Appendix we give the explicit expressions for 
the photon scattering rates used in Sec. IV following the 
formalism given in [61]. The total rate of photon scat¬ 
tering by an atom in the clock state ]/, 0) is 

7/=-^Im[(/,0l/ieff|/,0)], (Bl) 
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where the effective non-Hermitian light-shift Hamiltonian 
for one atom is 


We seek the time evolution of the one and two-body cor¬ 
relation functions: 


= (B2) 

as follows from Eq. (38), where ao(^//') = 

~McoA — i^ir /2 complex polarizability and 

the irreducible tensor operator A(/, /') is given in 
Eq. (46). 

The rate of optical pumping between clock states 
|/,0)^ |/,0) is 


Q 

where 27774^2 

Lindblad jump operators for optical pumping between 
ground levels f ^ f [61]. Each jump operator is 

associated with absorption of the probe photon polarized 
along ei„ followed by spontaneous emission of a photon 
with polarization e^, where q = {0,±1} labels spherical 
basis elements for tt and a± transitions. 

To find the dependence on the input field inten¬ 
sity, we define a characteristic photon scattering rate, 
7 s = with Rabi frequency H = 

■^3 *" 

2{j\\d\\j')£[^'^/h, reduced optical dipole matrix ele¬ 
ment (jjldjlj')) field amplitude = |E[j]^^(r')|. 

Eqs. (Bl) and (B3) yield. 


(iVc) = ^(ig‘)) 

n 

(Gla) 


(Clb) 

W = ^E(^ 3 ”^) 

(Glc) 

/ 92 \ (-^C) 1 1 /-(m) ..(n), 

(•^3) = 4 + 4 E ). 

(Cld) 


m^n 


where Ic = |t)(tl + IDdl is the single-atom projector 
onto the clock states. To include optical pumping, we 
apply the following equations of motion. For a collec¬ 
tive, single-body operator, X = evolution 

due to optical pumping is d(X)|op = Tr[I?„[y6]X]dt = 
[a;("i])c?f, where the map, which acts locally on 
atoms along the nanofiber, is given in Eq. (72). Two- 
body microscopic operators decay by optical pumping 
according to [i5] 



op 




(C2) 


where the superscripts refer to the and atoms. 

Applying the adjoint map to the single-atom operators 
yields 


7/ = ngN ]^ct ( A///) u 7„ (r() • (/,0 1 a (/, /' ) I /, 0) • Ui„ (r() 
/' 




/' 




E 






eg-{eine*n)-eg, 


(B4a) 

(B4b) 



e*.(ei„erj-e„ (B5) 


where a{Affi) = aoT^/AA'jj:, is the scattering cross sec¬ 
tion at the probe detuning, = {f'q \ /O; Iq) are the 

Clebsch-Gordan coefficients, and 

|oj/T = (2/ + 1)(2/ + 2)| I (B6) 


are the relative oscillator strengths determined by the 
relevant Wigner 6-J symbol. 


Appendix C: Derivation of the equations of motion 
for the moments 

In this Appendix we derive the equations of motion 
for the correlation functions that define the metrologi- 
cally relevant squeezing parameter, = NaAJ^/{Ji)^. 



[ic] = —7ooi c + 7 o3<5’3 

(G3a) 


= -7330'3 + 73oic 

(C3b) 


= -7iiO'i, 

(G3c) 

with rates 


700 = 

7T + 74 - 7T^T “ 7T^4 “ 74^4 “ 74^4 

(G4a) 

2 

703 = 

-7t + 74 + 7T-*-t + 7t-*-4 “ 74-*-4 “ 74-*-T 

(C4b) 

2 

733 = 

74 + 74 - 74^4 + 74^4 “ 74^4 + 74^4 

(G4c) 

2 

730 = 

-74 + 74 + 74^4 - 74^4 “ 74^4 + 74^4 

(C4d) 

2 

711 = 

74 + 74 

2 

(G4e) 


Given Eqs. (G2, G3), the equations for the two- 
body spin correlations, Eq. (77), follow. Similarly, one 
can derive equations of motion for the remaining two- 
body microscopic operator correlations G and 
(i^^ G -I- (Tg™ G ip when m ^ n and from these, 
the macroscopic operator expectation values (J|), (Nq), 
and (AcJg). As we have examined numerically, on the 
time scale of the QND measurement, the correlation be¬ 
tween atom number in the clock state subspace and the 
pseudospin moment is weak, and one can thus treat the 
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atom number operator in the clock state subspace as a 
c-number. We therefore set (iVc J3) — (-/Vc')( Js) = 0 and 
(N^) - (Nc)'^ = 0 and dehne Nc = {Nc)- 
The equations of motion for the moments of J 3 are now 
found from the SME, Eq. ( 68 ), 

d(J3) =sv^AJ| dW - 'y33{J3)dt + ^j3oNcdt, (C5a) 
d(j|) =2sv^( J 3 ) A J| dW -2733 {Ji)dt+l (2733 - 700 )Ncdt 
+ 'y3o{J3)Ncdt + ^ (703 — 2730 ) {J 3 )dt. (C5b) 


The stochastic term in d{J^) was simplihed by assuming 
Gaussian statistics [77], (Jf) = 3(J|)(J3) — 2{J3)^. Fi¬ 
nally, the Ito calculus governing the stochastically evolv¬ 
ing moments requires that differentials be taken to sec¬ 
ond order, and the evolution of the variance is given by 
dAJ| = d(J|) — 2{J3)d{J3) — (^(Js))^. The equation 
of motion for the conditional variance, Eq. (73d), then 
follows from Eqs. (C5). 


